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Abstract 

We propose a theory of low temperature thermal transport in nano-wires in the regime where a 
competition between phonon and flexural modes governs the relaxation processes. Starting with the 
standard kinetic equations for two different types of quasiparticles we derive a general expression 
for the coefficient of thermal conductivity. The underlying physics of thermal conductance is 
completely determined by the corresponding relaxation times, which can be calculated directly for 
any dispersion of quasiparticles depending on the size of a system. We show that if the considered 
relaxation mechanism is dominant, then at small wire diameters the temperature dependence of 
thermal conductivity experiences a crossover from T 1 / 2 to Independence. Quantitative analysis 
shows reasonable agreement with resent experimental results. 

PACS numbers: 63.20.Kr 63.22. +m 65.40.-b 
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Low-dimensional materials have attracted considerable attention in recent years, particu- 
larly, in view of their potential applications in electronic devices 0] . Many theoretical and 
experimental studies of nanowires and nanotubes are centered on the properties of electronic 
transport. However it is realized now that the thermal properties of nanomaterials are also 
important for applications 0,0,51. It is of special interest to increase thermal conductance 
in the micro and nanodevices jj, 6,(7]. In this paper study thermal transport in nonmetallic 
systems, in which heat is transported by thermal excitations only. In addition to practical 
importance of such studies, thermal transport in nanowires is interesting from fundamental 
point of view. Recent theoretic, fl and experiment^ findings proved the existence of 
the quantum of thermal conductance in ballistic regime, which is similar to the quantum of 
electronic conductance. The state of experimental and theoretical understanding of thermal 
transport in nanoscale systems is comprehensively discussed in the review Recently 



D. Li et al [ll| reported an accurate measurement of lattice thermal conductivity in silicon 
nanowires for a wide range of temperatures and wire diameters. They demonstrated the 
significant influence of the system size not only on the magnitude of the thermal conduc- 
tivity coefficient, but also on its temperature dependence. It is well known that for large 
enough diameters of the wire and diffusive phonon - boundary scattering, thermal conduc- 
tivity coefficient at low temperatures is proportional to T 3 . But for small values of the wire 



diameters experiment 11 1 shows clear crossover from cubic to near linear dependence on the 
temperature. In the present paper we consider one particular relaxation mechanism which 
can explain the observed crossover. 

Recen t l y Mingo Q carried out an accurate numerical s tud y of ther.ua! conductance of 
silicon nanowires to explain the decrease of the thermal conductivity coefficient with wire 
diameter observed in the experiment. He assumed that all the effects can be explained by the 
reconstruction of the phonon dispersion, where realistic phonon modes obtained from MD 
simulations were applied to general expression of the thermal conductivity coefficient. His 
numerical analysis shows excellent quantitative agreement with the experiment for large 
enough diameters at high temperatures. As the system size becomes smaller, the approach 
fails to describe a sharp decrease of thermal conductance as well as qualitative change of 
its temperature dependence. This is likely because the Matheissen's rule has been used for 
evaluation of the phonon lifetime, which has rather restricted range of applicability (see e.g. 
and references therein). 
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As it was noted in Ref. Q|, decrease of the temperature increases the characteristic 
phonon wavelength and reduces the scattering probability at the boundary surface. This 
leads to a modification of the phonon spectrum. For ideal wires it is represented by a 
set of branches with energies proportional to ID momentum directed along the wire. So 
the standard theory of thermal conductance in dielectrics and semiconductors has to be 
modified to account for low dimensionality effects as well as phonon spectrum modification 
at low temperatures. Thus, to understand thoroughly the physical processes occurring inside 
the nanowires with decreasing sizes, we need an analytical theory to account for different 
mechanisms explicitly, such as dispersion reconstruction and restricted geometry. 

To approach the problem we consider low enough temperatures, where the quasiparti- 
cle states of "acoustic" branches are thermally populated (e — > when p- — > 0). The 
corresponding acoustic branches have the following dispersion relations jl^ll"!^: 

ei = uxpi, e 2 = u 2 apl, (1) 

where q stands for the energy of a quasiparticle, pi is the corresponding momentum, a is 
the wire diameter, u\ and u 2 are the characteristic velocities. The first expression in Eq. 
is the phonon dispersion and the second expression is the dispersion of flexural mode. The 
nature of flexural modes comes from the fact they are analogous to bending modes of classical 
elasticity theory, or the antisymmetric Lamb waves of a free plate [16] . Appearance of such a 
mode is just a direct consequence of restricted geometry, and under some conditions it can be 
considered as the only size effect in thermal transport properties. Strictly speaking there are 
two modes for each type of dispersion, but we do not account for them separately since their 
contributions are qualitatively the same. Consequently, we have to solve the kinetic problem 
for two-component gas of quasiparticles. It is well known jl^, 12, 3,13,0] that using simple 



Callaway formula to estimate the thermal conductivity coefficient in two-component systems 
sometimes leads to serious confusions. Simple summation of the relaxation rates, as is done 
in the majority of theoretical works, is questionable under many physical conditions. The 
dependence of the kinetic coefficients on different relaxation times is much more complicated 
in reality. Accurate method for calculation of the diffusion coefficient in two-component gas 
of quasiparticles was proposed in Ref. [^j. Here we extend this formalism to thermal 
conductance problem. 

To start, we consider a system of two types of quasiparticles. Their kinetics is described 



by equations for corresponding distribution functions ff. 

dz 



vM = J2c ij (f i J j ) + C i3 (f i ), 2 = 1,2. (2) 



where Cjj(/j, fj) is the collision integral of thermal excitations, and Cis(fi) is the collision 
integral describing the scattering processes between quasiparticles and scatterers. In gen- 
eral the latter includes all the processes leading to the non-conservation of the total energy. 
Vi = dei/dpi is the group velocity of the corresponding thermal excitation. The main purpose 
of our theory is to obtain analytic expressions of thermal conductance, which are applicable 
to quasiparticles with arbitrary dispersion relations. In other words, the explicit dispersion 
relations in Eq. are needed only at the last stage when calculating corresponding relax- 
ation times and thermodynamic quantities. As usual, we seek a perturbative solution of the 
system (2) in the form 

fi = if + Sf it (3) 
where is the local equilibrium Bose-function and 5fi <C represents small devi- 
ation from the equilibrium. The perturbation term can be conveniently chosen to be 
5fi = —gidfi /dei with g^ the new target functions. After the standard linearization pro- 
cedure, Eq. (J2J) can be written in the following matrix form: 

1 dT 

l**>r^- = <%>. ( 4 ) 

where 

eivi 



\g) 



The 2D collision matrix C can be decomposed into a sum of three terms, corresponding 
to different relaxation mechanisms: C = J + S + U, where J" , with matrix elements Jij = 
CikSij + Cij(l — 5ij) (k 7^ i), describes the relaxation due to interaction between quasiparticles 
of different types; S = CjAj) describes collisions between identical quasiparticles, and 
U (Uij = Ci^5ij) describes all the other relaxation mechanisms, which do not conserve total 
energy of the quasiparticle system, hi includes scattering on defects, boundaries, umklapp 
processes etc. |2l( Here Cy represent linearized collision operators. 

Let us define the scalar product of two-dimensional bra- and ket-vectors as follows 

(01%) = E MX*) = - E / <PlXk^dT k , (5) 



where (<fik\ and \xu) are the correspondent one- component vectors, dT is the element of 
phase volume. Under this condition, the collision operator C becomes hermitian. System 
(J3J) is the system of nonuniform linear integral equations. According to the general theory of 
integral equations the target solution \g) must be orthogonal to the solution of corresponding 
uniform equations C\(f) un i). It is therefore convenient to write the formal solution of (@J) so 
that the orthogonality condition (g\4> un i) is imposed explicitly in the solution. For this 
purpose we define the projection operator T n onto the subspace orthogonal to the vector 
\4>uni), V n = 1 — Vet V c = \4>uni){4>uni\- As a result, the formal solution of the system (JJJ) 
can be written in the form 

\ 9 } = A{c- 1 )vM^. (6) 

The heat flux density due to the thermal excitations of different types is given by the 
expression Q = J2k=i 2 / e kVkfkdT k . Using relation (jSJ) and definition of scalar product 
Q can be rewritten as Q = ((j) K \g). On the other hand, the effective thermal conductivity 
coefficient is defined by the relation Q = —K e ffdT/dz. Comparing the above two expressions 
for Q, and using the formal solution (JUJ) we obtain 
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ye- 1 !^}. (7) 

To derive an exact and analytical expression for thermal conductivity coefficient (J7J) it is 
necessary to introduce a complete set of orthonormal two-dimensional vectors \ip n ) (n = 
1,2,3,...) belonging to the infinite-dimensional Hilbert space with scalar product (JSJ). In 
principle, the particular choice of the basis is not essential, but for the convenience of 
calculations it is useful to specify at least four of them. It's convenient to chose the first 
of them to correspond to the total momentum of quasiparticles and the second one to be 
orthogonal, but still linear in momentum 



iv*> = 4= 




IV*> 




(8) 



\fPP1P2 

where pi = (pi\pi) is the normal density of the ith component, p = p\ + p^- The third and 
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the fourth vectors correspond to the energy flux: 



where 




1^4) 



t<2 




(9) 



If ST 

= VT \ jVj ~ ~fr Pj ' ' ^ 



and M K j = \j (il>Kj\4>nj) is the corresponding normalization coefficient. Partial entropy of 
quasiparticle subsystem Sj in eq. (|1U|) is given by the relation 

S J = f( e 3 V j\Pj)- (11) 

Formally, the kinetic problem of two-component quasiparticles system can be solved in 
the above basis set. The inversion of the operator matrix C in Eq. (171) is similar to the 

n 

procedure described in Ref . |20J . The final result contains infinite dimensional non-diagonal 
matrices. To obtain closed form expressions we must use some approximations, correct r- 
approximation [19] or Kihara approximation 0, 0, Q- 

In some physical situations we 
are able to obtain closed analytical expressions. It is rigorously proved in that in case 
of quasi-equilibrium within each subsystem of quasiparticles, the corresponding transport 
coefficient can be obtained in close analytical form. This is a reliable approximation when 
the low temperature relaxation is mainly governed by the defect scattering processes. The 
approximation formally implies that all the matrix elements of matrix S in eq. tend to 
infinity. The thermal conductivity coefficient in this case can be obtained in the form: K e ff = 
Kp + kd- Here we separate the flux part of thermal conductivity coefficient Kp = tfS 2 T / p 
with S = Si + S 2 , which approaches infinity when the quasiparticles do not interact with 
scatterers, and the diffusive part k d = t d (S\T / pi — S 2 T / p 2 ) 2 p\P2/Tp. The corresponding 
relaxation times are given by 

' D = \ 

and 



t» = ^t£ + Jt£+t£ + t£} , <12) 
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Relaxation times contained in formulas (|12|) - are defined by 



V = ^(Pil^bi)- (i4) 



We emphasize that they are not actual scattering times, which are momentum dependent, 
but are relaxation times associated with corresponding scattering mechanisms. Once we 
obtain the particular scattering rate fkj{Pk) from standard scattering theory, we can replace 
true collision operator C k with Vkjipk), so that the corresponding relaxation time can be 
calculated by 



(0) 



kj - Pk J Pk^j(Pk)^-dT k (15) 

As can be seen from the derived formulas, the coefficient of thermal conductivity contains 
different relaxation times in rather non-trivial combination. If drop out one component (say 
set S 2 = 0, P2 = 0) we recover the usual result = t V3 S 2 T / pi. For phonons with linear 
dispersion e = vp, k f reduces to the well-known result K P p = C p / l f 2 r 13 /3, where C p h = 3S p h 
is the heat capacity of phonon gas. 

The main advantage of our approach is its universality. In fact till now we have not 
restricted ourselves to any particular dimensionality of the system or any quasiparticles 
dispersion. All the necessary information is contained in the corresponding relaxation times 
and thermodynamic quantities. This formalism allows us to analyze contributions from 
different relaxation mechanisms to the total thermal conductivity coefficient. Given the 
dispersion relations of quasiparticles, we can easily calculate all the quantities contained in 

With Eqs. (fT2j) - (fT3Jl we are able to address the competition between relaxation processes 
of the flexural and phonon modes. Glavin 14] noted that such a competition can be essential 
at extremely low temperatures if the dominant relaxation mechanism is elastic scattering 
on defects, where he argued that thermal conductivity coefficient would scale as T 1 / 2 . Our 
approach allows us to study this competition comprehensively. Particularly, we predict the 
strong dependence of the temperature scaling exponent on the wire diameter. The standard 
Fermi golden rule approach [l4| gives the momentum dependent scattering rates for different 

1/2 

modes v 13 = ^13:^372, v 23 = ^23^3, where W kj are the corresponding scattering 
amplitudes, which depend on the physical properties of particular material. Using Eq. (fTKj) 
it is easy to show, that the corresponding relaxation times scale as 

oc a-^T 1 / 2 , r^ 1 oc a-^T- 1 / 2 . (16) 



Different temperature dependence of relaxation times lead to a strong competition be- 
tween two physically different mechanisms of thermal conductivity, flux and diffusive. The 
dominance of one over the other strongly depends on the wire diameter at a given temper- 
ature. To make some specific conclusions let us summarize the approximations done and 
specify the range of validity of the proposed theory. We consider a situation when thermal 
excitations are multiply scattered elastically while being transferred through the wire, so 
that other scattering mechanisms are strongly suppressed by interaction with defects. Only 
for this case we were able to drop relaxation within each subsystem of identical quasiparticles 
to obtain closed expressions (JEJ), (JXHJ) - The influence of boundary is accounted in the dis- 
persion of flexural mode and in the dimensionality of the system. The range of temperature 
is supposed to satisfy the relation T < Ae, where Ae ~ 1/a is the characteristic value of the 
frequency gap between the adjacent phonon branches. For larger temperature we cannot 
use the acoustic modes (HJ only, but need to account higher branches. 

In Figure 1 we compare our theoretical results with the experimental data from |TJj . We 
have chosen the unknown parameters W13 = 1.2 • 10 _44 m 5 s -4 and W23 = 0.9 • 10 _44 m 5 s -4 
to fit data for a = 22nm. Deviations from the experimental data for large diameters and 
temperatures show the restriction of applicability of our initial approximations. They arise 
from the Debye approximation and simplified dispersion expression. Additionally, when the 
diameter of the wire increases, the mechanism we consider becomes less dominant. To be 
more precise we need to include higher excitation branches as well as another relaxation 
mechanisms. However our approach allows to understand the physics of the processes in the 
region under consideration. Clearly, observed crossover is the result of competition between 
Kp oc T 1 / 2 and Kp, oc T 3 . For smaller diameters, Kp is strongly dominant in the wide range 
of temperatures as shown on Figure 2 for wire diameter a = 2nm. Figure 3 demonstrates 
a complete crossover from T 1//2 to T 3 dependence for the nanowire of a = 30nm. It can be 
seen that the T dependence between 20K and 40K is nearly linear, which was observed in 
the experiment It should be noted that T 3 dependence of Kp> cannot be interpreted 

in simple analogy with the bulk case. It comes not from a specific heat directly, but from 
different sources including competition of the relaxation times in Eqs. (II 2|) . ()13|) . 

In summary, we derived the general analytical expressions ()12|) - (fT3*j) . to explicitly cal- 
culate the contributions of different scattering mechanisms to the total relaxation of the 
system. The simple expressions clarify the essential effects leading to the observed behavior 
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of thermal conductivity coefficient. It is clear that the particular dispersion laws (and their 
reconstruction) affect scattering rates and thermodynamic quantities. Restricted geome- 
try and low dimensionality lead to additional scattering mechanisms. Note the information 
about dimensionality is naturally included in the particular form of phase space element oTV 
Such a formalism helps to distinguish effects from different scattering mechanisms. When 
applied to the regime, where phonon modes compete with flexural ones, our theory agrees 
favorably with available experimental data. Furthermore, we showed that the thermal con- 
ductivity coefficient changes from approximately T 1 / 2 -dependence to Independence with 
increasing temperature. In view of our theoretical results it is useful to investigate smaller 
diameters or lower temperatures with fixed diameters in experiment to better reveal the 
crossover from T 3 to T 1 / 2 dependence. 

We would like to thank Prof. Peidong Yang and Dr. Deyu Li for providing us the 
experimental data. This research is supported by the NSF Career Award (Grant No. Che- 
0093210). J. C. is a recipient of the Camille Dreyfus Teacher-Scholar Award. 
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FIG. 1: Thermal conductivity coefficient calculated from Eqs. (|12|) and (jlHj) for different values of 
nanowire diameter. Experimental data are from Ref. 
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FIG. 2: Comparative contribution from flux and diffusive parts of thermal conductivity for a 2nm 



wire. 
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FIG. 3: Crossover from T 1 / 2 to T 3 dependence of K e // for a 30nm wire. 
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